Transient tunneling effects of resonance doublets in triple barrier systems 
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Transient tunneling effects in triple barrier systems are investigated by considering a time- 
dependent solution to the Schrodinger equation with a cutoff wave initial condition. We derive 
a two-level formula for incidence energies E near the first resonance doublet of the system. Based 
on that expression we find that the probability density along the internal region of the potential, is 
governed by three oscillation frequencies: one of them refers to the well known Bohr frequency, given 
in terms of the first and second resonance energies of the doublet, and the two others, represent a 
coupling with the incidence energy E. This allows to manipulate the above frequencies to control 
the tunneling transient behavior of the probability density in the short-time regime. 

PACS numbers; PACS; 03.65. Xp,73.40.Gk 



I. INTRODUCTION 



In this work we address the issue of time-dependent 
tunneling phenomena in triple barrier resonant systems, 
aimed to study the transient behavior of the probabil- 
ity density near a resonance doublet. We shall refer to 
these structures as two-level open systems, in the sense 
that their finite- width barriers enable the system to inter- 
act with incident particles via a tunneling process. The 
dynamical properties of triple barrier structures have not 
drawn the attention they deserve. In this work we wish to 
emphasize that triple barriers involve novel dynamical as- 
pects not present in double-barrier structures, where the 
tunneling dynamicSunaax resonance energy is governed by 
a single resonance. Elalj^a 

The purpose of this paper is to demonstrate, based on 
an exact analytical approach, that the dynamics of the 
transient probability density is governed by three rele- 
vant frequencies that involve the resonance energies of 
the doublet and the incidence energy E. We find that in 
addition to the Bohr frequency, LO21 = \£i — ^il/K which 
is an intrinsic property of the system, there are two addi- 
tional frequencies, ui = \E — £i\/h and uj2 = \E — £2\/fi, 
where the resonance energies £1 and £2 are the real 
parts of the corresponding complex resonance energies 
En = £n — i^nf^ (n=l,2) of the problem. This should 
be contrasted with the well knowntl dynamical behavior 
of closed two-level systems, which is governed only by 
the Bohr frequency, fJi2 = (E2 — Ei)/S where Ei and 
E2 are the real energy eigenvalues of the system. As 
shown below, the above frequencies may be manipulated 
to produce a significant enhancement of the short-time 
transient behavior of the probability density. 

The paper is organized as follows: Section II presents 
an overview of the formalism. In section III we discuss 
the transient behavior of the probability density through 
several numerical examples. Finally, section IV, provides 
some concluding remarks. 



II. THE FORMALISM 

The model used in this work deals with an explicit so- 
lution to the time-depeudent Schrodinger equation with 
cutoff initial conditions ,d and consists of a generalization 
to tunneling problems of the free auantum shutter setup 
that predicts diffraction in time!Si!M The phenomenon 
of diffraction of matter in time has been recently ex- 
perimentally verifiedBEj and has also stimulated further 
studies .Ell The setup used in this work may be visualized 
as a quantum shutteiil3 placed at a; = 0, just at the left 
edge of the resonant structure that extends e«er the in- 
terval < X < L. Upon opening the shutterll3 at t = 0, 
the incoming initial wave, represented by a cutoff plane 
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x > 0, 
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interacts with the internal region (0 < a; < L) of the po- 
tential. The wave solution to the time-dependent prob- 
lem "^{x, k; t) for a; > and t > 0, is given by,EI 



* = $feAf(yfc)-$_feM(y_fe) 



E 



PnMivkJ. (2) 



The quantities ^±k = ^{x, ztk) refer to the stationary 
wave solution, and the factors. 



Pn{x, k) = 2ikUn{Q)Un{x) / (fc^ - fc^). 



(3) 



are given in terms of the resonant states {un{x)} and the 
complex energy eigenvalues E^ — h?kn/'2m of the prob- 
lem. The complex energy eigenvalues may be written 
in terms of the complex wave numbers kn — ttn — ^^n: 
and correspond to the S-matrix poles of the problem. 
They are distributed in the third and fourth quadrants 
on the complex fc-plane iuj-s, well known manner. The 
A'/-functions are defined as,c 



1 



(4) 



2 



where the-.functions w{iys) stand for the complex error 
functionJii 



w{iys) = e^=eTk{ys) 
and the argument ys, reads, 

,..e-/M_i.]^^,V. 



(5) 



(6) 



and s stands for ±fc or fc±„. As shown elsewherep the 
time-dependent solution given by Eq. ^ goes into the 
stationary solution $fc at asymptotically long times. 

For triple barrier systems, the resonance spectra typi- 
cally corresponds to a succession of resonance doublets, 
formed by the coupling of the single resonances associ- 
ated with each of the two wells of the system. We shall be 
interested in systems where the first doublet is isolated. 
The approximation to Eq. (||) then reads, 

2 

-Y,{pnM{ykJ+P-nM{yk_J} , (7) 



where we have used $_fc — <i>^. For a resonance doublet 
the stationary function may also be. written as the sum 
over the first two resonance terms, lj namely. 



$fc(x, k) « pi{x, k) + p2(x, k), 



and consequently. 



\Pi\ 



\P2\ 



Pl2, 



(8) 



(9) 



where pi2 = 2Ke{pip2}. Although the time-dependence 
of Eq. (|^) is contained in the i\f -functions, a considerable 
simplification of this two-level formula can be derived, in 
which the time dependence is explicitly given in terms of 
simple fuBGtions. Such a derivation is discussed in detail 
elsewherell3 and we will recount it here briefly. The M- 
functions M{yk) and M{yk^) contained in Eq. (0), can 
be related to functions of the form M(?/_fc) and M{yk_„) 
by means of the symmetry relation,B 



Using Eq. (|lo| ) we can rewrite Eq. (Q) as. 



(10) 



+ Aix,t), (11) 



where A{x,t) accounts for all the terms containing M 
functions of the form M{y^k) and M{yk_„), which be- 
have as aruinverse power of t, as follows from its series 
expansion,! M{y,) - l/2[l/(7ri/2y,) - l/in^/^y^J + ...]. 
Thus, except for extremely short or very long times com- 
pared with the lifetimes of the resonance levels of the 
doublet, the term A(x,t) gives a negligible contribution 
to the solution and can be neglected. By doing this, we 



can obtain a simple expression for the probability den- 
sity, valid for the internal region and an energy E close 
to the doublet, namely. 



\^iE,t)f = ME,t) + ME,t) + ^i2{E,t), (12) 

where (j)n{E,t) and the interference terms (j)mn{E,t) 
(n=l,2) are given respectively by. 



(j)n{E,t) = \pn\ Xn{E,t), 



(13) 



and 



^^n{E,t) = 2Re{pmP:,Un{E,t)}, (14) 



where E — h'^k^/2m is the incidence energy, and the 
functions Xn and ^„i„ have the following closed analytic 
expressions, 

XniE, t) = l-2 cos{uJnt)e-^-"^'' + e-^"*/^ (15) 



(16) 



In the above expressions, (Di, UJ2 and 0)12 are defined by 
LJn ^ {E - £n)/fi and 0)21 = (£2 - Si)/fi- 

The formula given by Eq. ( [T^ ) is an important an- 
alytical result since it explicitly reveals novel aspects 
of the quantum dynamics of tunneling structures with 
resonance doublets. According to Eq. (^2|), the time- 
dependent probability density is the superposition of 
the three oscillating contributions, (j)i{E,t)^ 4)2{E,t) and 
(j)2i{E, t), which have in general different amplitudes and 
frequencies. The three characteristic frequencies that 
govern the time evolution during the transient regime 
are: uji, uj2 and 0^21, which are respectively the absolute 
values of tDi, 0)2 and uj2i- Note that at asymptotically 
long times, it is easily seen from Eqs. (15) and ([l^), re- 
spectively, that Xn ^ 1 and ^mn 1, and hence the 
probability density for the two-level formula, given by 
Eq. (12), goes into the stationary solution given by Eq. 

(I)- 



III. EXAMPLES 

We shall be interested in analyzing the transient tun- 
neling effects of the probability density at the right-hand 
edge of the system, x = L, because there appear the 
largest transient effects along the transmitted region. We 
consider as a first example a periodic triple barrier sys- 
tem with parameters given as in Ref. |l^, namely: bar- 
rier heights Vo = 0.12 eV, barrier widths 60 = 3.0 nm, 
well widths wq — 16.0 nm; and effective mass of the 
electron m — 0.067me. The corresponding resonance 
parameters of the first doublet are: energy positions, 
£1 = 11.512 meV and £2 — 14.387 meV; and resonance 
widths, Fi = 0.4089 meV and r2 = 0.6365 meV. 
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FIG. 1: Time evolution of \'i'{L,k;t)f for a triple barrier 
system at off-incidence energy E — £'i-|-2.0ri, using the exact 
solution, Eq. M with A*' = 4 (solid line), and the two level 
formula, Eq. ([L2|) (dashed line). The systems parameters are 
given in the text. 



Let us first illustrate the reliability of our approximate 
formula derived above for a single doublet. In Fig. |l| 
we compare the behavior of the probability density using 
both the formal solution, Eq. (^) (solid line), and Eq. 
( |l^ ) (dashed line), for an incidence energy near the first 
resonance, E = £i + 2.0Ti = 12.33 meV. As can be ap- 
preciated, the two level approximation to Eq. (H) given 
by Eq. (|l2|) gives an excellent description of the probabil- 
ity density. In this particular example, we have included 
in Eq. (J^ the first four resonances of the systems, i.e. 

= 4, in order to illustrate that the contribution of 
far away resonances is negligible. The irregular behav- 
ior of l^tp observed in Fig. arises from the interplay 
between (pi, 02, and (j)i2 of Eq. ([l^). This situation 
contrasts with the regular behavior observed iii. double 
barrier structures. As shown in a recent work,a in the 
case of a double barrier system the probability density 
grows exponentially for incidence at resonance, and ex- 
hibits regular oscillations with a. single frequency if the 
incidence occurs near resonance.EI This is due to the fact 
that in the double barrier case, the one-level approxima- 
tion stands, and hence only the term (j)i{E,t) is impor- 
tant. An interesting feature of triple barrier systems, not 
present in double barrier structures, is that the frequen- 
cies can be tuned by a proper choice of the incidence 
energy E. This allows to manipulate the frequencies in 
such a way that the irregularities observed in Fig. |l| diss- 
appear. This occurs at two special situations that depend 
on E. The first situation is when the incidence energy 
coincides with one of the two resonances, and the second 
one occurs when the incidence energy coincides with the 
middle point between the two resonances of the doublet. 
In the first case, only one of the three terms of Eq. ( |l^ ) 
dominates over the remaining two, for example li E — £i, 




30.0 



FIG. 2: Time evolution of [^(L, fc; t) |^ for the same triple 
barrier system of the previous figure, for two special situ- 
ations: (a) for incidence at the first resonance, E = £i, 
where a>i = and a>2 ~ —1^21; and (b) for incidence at 
E = E = {£1 + £2)72, where tii = —lu2 = l/2ii2i. For com- 
parison, the calculations in (a) were made by Eq. (|7|) (solid 
line) and the exponential formula given by Eq. (|l7|) (dashed 
line); and the calculations in (b), by Eq. ^ (solid line) and 
Eq. (Ili) (dashed line). 



then (j)2{E,t) and 4'2i{E,t) are negligible in comparison 
with (f)i{E,t). Since in this case lui = 0, then wc have 
that the probability density is governed by the following 
simple expression. 



|*(^-£i)|2 «T(£i)(l 



-t/ri^2 



(17) 



where ri = h/T is the lifetime of the resonance n — 1 and 
T{£i) is the peak value of the transmission coefficient, 
which is unity for this symmetrical system. The results 
of this resonant case are depicted in Fig. |^(a), where 
we show the calculations using Equations (0) (solid line) 
and (1^) (dashed line). Both curves almost coincide, ex- 
cept by the very small oscillations of the exact two-level 
formula, i.e. Eq. (^ (solid line), which are due to the 
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FIG. 3: (a) Enhancement of k; t)p (solid-line) in a triple 

barrier systems for incidence at the center of the resonance 
doublet, compared with the results of a double barrier system 
(dashed- line), (b) Enhancement as a function of the central 
barrier width b = 3.0 nm (solid line), b = 4.0 nm (dashed- 
line), and b = 5.0 nm (dotted- line). The arrows indicate 
for each case the values of the corresponding transmission 
coefficient. See text. 



width wq = 5.0 nm. The first resonant state of the sys- 
tem has energy position, £i = 80.11 meV; and resonance 
width, Fi — 1.033 meV. The incidence energy was also 
chosen with the same deviation from resonance in units 
of the resonance width, that is, E = £i + 3.515ri whose 
numerical value is 83.740 meV. Note that the scale in the 
time axis was normalized to the lifetime of the first reso- 
nance of each system, which for the triple barrier has the 
value Ti = 1.61 ps, and for the double barrier, n = 6.37 
ps. Both curves tends to their correct asymptotic limits, 
the transmission coefHcient, which for the double barrier 
system has the value T(E)=0.0229, and for the triple bar- 
rier system, T(E)=0.119. These values are indicated by 
the small arrows in Fig. ||(a). 

As it is well known p^cap time-independent studies in 
triple barrier systems Jl3ll3 for incidence energies at the 
center of an isolated doublet, the transmission coefficient 
increases to unity as we increase the width &2 of the cen- 
tral barrier to about twice the value of the width of the 
external barriers. In view that the transmission coeffi- 
cient is the asymptotic value of the time dependent prob- 
ability density at x — L, it is expected that the latter 
can also be enhanced in the same fashion. In Fig. ^(b) 
we illustrate how we can manipulate this extra degree 
of freedom to enhance the amplitude of the oscillations 
of the transient probability density. Here we considered 
&2 = 4.0 nm (dashed line), and &2 — 5.0 nm (dotted line); 
the solid line corresponds to the same triple barrier sys- 
tem of Fig. |(a) (62 = 3.0 nm), which is also included 
here for comparison. Also in this figure, the values of the 
transmission coefficient at the energies E for each system 
is indicated by arrows to illustrate how the time depen- 
dent probability density tends to the correct asymptotic 
behavior as t —>■ 00. 



IV. CONCLUDING REMARKS 



effect of the second resonance of the doublet. In the sec- 
ond case, when the incidence energy is chosen just at the 
middle of the two resonances i.e., E = E = {£1 + f2)/2, 
we have w = cDi = —0)2 = 1^21/2, that is, the dynamics 
is governed b;^ a single frequency, (D, and the behavior 
of I'i' (Zra/c; t)| vs t is similar to a diffraction in time 
pattern,Ei see Fig. §(b). Here the numerical value of 
E = 12.949 meV= £1 + 3.515ri. 

Compared to the double barrier case, this constructive 
interference effect produces an important enhancement 
of the transient probability density. The comparison is 
shown in Fig. ^(a) in which we used the same triple bar- 
rier structure parameters of the previous figures, and the 
double barrier system with potential parameters: barrier 
heights Vq = 0.23 eV, barrier widths bo = 5.0 nm, well 



In conclusion, the dynamics of the probability density 
for triple barrier resonant structures, which is a typi- 
cal example of an open two-level system, has been ex- 
plored. We have derived a simple analytic expression 
for the probability density that provides an accurate de- 
scription for energies near the resonance doublet of the 
system. The two-level formula allows to identify three 
relevant frequencies that govern the transient behavior 
as a function of time. The derived formula goes into 
the stationary two-level solution at asymptotically long 
times, and thus establishes a link with the usual sta- 
tionary approach. Our results suggest that the transient 
effects that we have discussed are of relevance at short 
times and distances from the interaction region. Hope- 
fully our results may stimulate experimental work on this 
subject. 
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